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Abstract 

Including contributions of scale-dependent vacuum expectation values, we 
derive new analytic formulas and obtain substantially different numerical pre- 
dictions for the running masses of quarks and charged-leptons at higher scales 
in the SM, 2HDM and MSSM. These formulas exhibit significantly different 
behaviours with respect to their dependence on gauge and Yukawa couplings 
than those derived earlier. At one-loop level the masses of the first two genera- 
tions are found to be independent of Yukawa couplings of the third generation 
in all the three effective theories in the small mixing limit. Analytic formulas 
are also obtained for running tan in 2HDM and MSSM. Other numer- 
ical analyses include study of the third generation masses at high scales as 
functions of low-energy values of tan /3 and SUSY scale Ms = Mz — 10^ GeV. 
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I. INTRODUCTION 



One of the most attractive features of current investigations in gauge theories is the 
remarkable unification of the gauge couphngs of the standard model (SM) at the SUSY 
GUT scale, Mjj = 2 x 10^^ GeV, when extrapolated through the minimal supersymmetric 
standard model (MSSM) [|T|. Although the nonsupersymmetric standard model (SM), or the 
two-Higgs doublet model (2HDM) do not answer the question of gauge hierarchy, unification 
of the gauge couplings is also possible at the corresponding GUT scales when they are 
embedded in nonSUSY theories like SO{10) and the symmetry breaking takes place in two 
steps with left-right models as intermediate gauge symmetries @. Grand unification of 
gauge couplings of the SM in single-step breakings of GUTs has also been observed when 
the grand desert contains additional scalar degrees of freedom and the minimal example 
is a ^(3, 0, 8) of SM contained in 75 C SU{5) or 21Q C 50(10) with mass = - lO^^ 
GeV [0. Unification of gauge couplings in nonSUSY 5*0(10) has been demonstrated with 
relatively large GUT-threshold effects 0. Yukawa coupling unification at the intermediate 
scale has also been observed in nonSUSY 50(10) with 2HDM as the weak scale effective 
gauge theory [Q. Apart from unity of forces at high scales, SM, 2HDM and MSSM have 
tremendous current importance as effective theories as they emerge from a large class of 
fundamental theories. 

Recent experimental evidences in favour of neutrino masses and mixings have triggered 
an outburst of models many of which require running masses and mixings of quarks and 
charged- leptons at high scales as inputs for obtaining predictions in the neutrino sector . 
The running masses are not only essential at the weak scale, but they are also required 
at the intermediate and the GUT scales in order to testify theories based upon quark- 
lepton unification with different Yukawa textures and for providing unified explanation of 
all fermion masses p|-p!^. Quite recently extrapolation of running masses and couplings to 
high scales have been emphasized as an essential requirement for testing more fundamental 



theories 12 
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In a recent paper one of us (M.K.P) and Purkayastha have obtained new analytic 
formulas and numerical estimations for the fermion masses at higher scales in MSSM includ- 
ing contributions of scale-dependent vacuum expectation values (VEVs) where the SUSY 
scale (Ms) was assumed to be close to the weak scale (M5 ^ Mz). In this paper we extend 
such investigations to SM, 2HDM and MSSM with the SUSY scale Ms > O (TeV). 

It is also possible that in a different renormalisation scheme, similar to that formulated 
by Sirlin et. al. fl^, the VEVs themselves do not run when they are expressed in terms 



of physical parameters defined on the mass shell. This makes it possible to avoid separate 
running of the VEVs and Yukawa couplings, but have just the fermion masses directly as 
running quantities. While it would be quite interesting to examine the consequences of such 



a scheme, the purpose of the present and the recent works [|T3| is to address the outcome of 
the most frequently exploited renormalisation scheme where the Yukawa couplings and the 
VEVs run separately flf^. 

This paper is organised in the following manner. In Sec. |1| we cite examples where 
running VEVs have been exploited by a number of authors and state relevant renormalisation 
group equations (RGEs). In Sec. |ITT| we derive analytic formulas. In Sec. |V| we show how 
the formulas derived earlier for MSSM are modified when Ms ^ Mz- Numerical predictions 
at higher scales are reported in Sec. 0. Summary and conclusions are stated in Sec. 

II. RGES FOR COUPLINGS AND VACUUM EXPECTATION VALUES 

After the pioneering discovery of 6 — r unification at the nonSUSY SU (5) GUT scale 



24 1, a number of theoretical investigations have been made to examine the behaviour of 



Yukawa couplings and running masses at higher scales. Following the frequently exploited 



renormalisation scheme [15-23] where the Yukawa couplings and the VEVs run separately. 



the running Dirac mass of the fermion 'a' is defined as, 

MM=YMvaifi) (2.1) 
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Then the running of Ma{^) is governed by the RGEs of Ya{fi) and fa(/i) both. To cite some 
examples, Grimus has derived approximate analytic formulas in SM for all values of n 
extending upto the nonSUSY SU (5) GUT scale utilising the corresponding scale-dependent 
VEV. In the discovery of fixed point of Yukawa couplings, Pendleton and Ross have 
exploited the RGE of the SM Higgs VEV to derive the RGEs of the running masses from 
/X = Mw — Mgut- Anomalous dimensions occurring in the RGEs of respective VEVs 
have been explicitly derived and stated up to two-loops by Arason et. al. [1T5| , [T6| and by 



Castano, Pirad and Ramond for SM and MSSM. While investigating renormalisation 



of the neutrino mass operator, Babu, Leung and Pantaleone [23| have derived RGE for 



tan l3{fi) in a class of 2HDM as a consequence of running VEVs in the model. More recently 



Balzeleit et. al. |19] have utilised the RGE of the VEV in SM to determine running masses 



for fi = Mw - 10^° GeV. Cvetic, Hwang and Kim |20| have derived RGEs for the VEVs 



in 2HDM and utilised them to obtain running quark-lepton masses at high scales and also 
investigate suppression of flavour changing neutral current in the model. Most recently the 
RGEs of running VEVs have been utilised by one of us (MKP) and Purkayastha |TB[ who 



have obtained new analytic formulas and numerical estimations of the fermion masses at 
higher scales taking the SUSY scale Ms ~ Mz- 

We consider only the class of 2HDM where gives masses to up-quarks and to down- 
quarks and charged-leptons. For the sake of simplicity we ignore neutrino mass in the present 
paper which will be addressed separately. Our definitions and conventions for the Yukawa 
couplings and masses are governed by the following Yukawa Lagrangian (Superpotential) in 
SM or 2HDM (MSSM) and the corresponding VEVs of Higgs scalars, 
SM 

jOy = QlYu^Uu + QlYd<^Dr + IlYe^Er + h.c. 

(<|.°(^)) = v{^,) (2.2) 

2HDM, MSSM 

Cy = QlYu^uUr + QLYo^dDR + iLYE^dER + h.c. 



i'^dif^)) = Mt^) = cos (3{fi) 



(2.3) 



The relevant RGEs for the Yukawa matrices at one-loop level for the three effective theories 



are expressed as [p^5Hl8| , |25| --P7 
2dYu 



167r 
167r 
167r 



dt 
^ dYp 
dt 

2 dYE 
IT 



Tr {sYuYi + SaYoYl, + ar^r^) + - {bYuY^ + cYdY^) - ^ C^^g^ 
Tr {SaYuY^ + SYpY^ + YeY^) + | {bYoY^ + c^/FJ) - J] C^'^g^ 



Tr {3aYuY^ + 3YnY^ + YeY^) + -BYeY^ - ^ ^^^^ 



Yr 



Yu 
Yd 

(2.4) 



The RGEs for the VEV in the SM has been derived up to two-loop from wave function 



renormalisation of the scalar field 0,|TB,|T^,|T^,pl|,E3] and the one-loop equation is 



,dv 



167r'^ = [E C^a' - Tr {sYuY^ + 3YnYl + YeY^ 



(2.5) 



where t = In /i. 

The RGEs for Va{a = u,d) in the 2HDM up to one-loop and in MSSM upto two-loops 
have been derived in fTE HTSl , ^ . The one-loop equations in both the theories are, 

,dvit 

'dt 
,dvd 



16n^^=[£C^g^-Tr{3YuY^ 



16^^ = Clgi - Tr (sFz.^^ + YeY^^ 



(2.6) 



The gauge couplings in the three models obey the well known one-loop RGEs, 

IGtt — = kgf 



dt 



(2.7) 



Two-loop contributions have been derived by a number of authors [Il5|- [lqj21| --p7| . The coef- 
ficients appearing in the RHS of (^^ )- (p.7| ) are defined in the three different cases, 
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SM, 2HDM 



MSSM 



SM 



2HDM 



MSSM 



20' 4' 
1 9 

4' I' 

Cr=(|;.^.0) (2.8) 
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16 




y. 


^■^^ 


16 


T, 


|3.0 


) 


3 3 




20' 4' 



ct- 

C:=[^r-,M (2.9) 



' VlO' 6 ' 

(a,6,c) =(1,1,-1) (2.10) 



6i= (y,-3,-7 
(a,6,c) = (o,l,^) (2.11) 



&.= (f,l,-3 
(a,6,c) = (o,2,^) (2.12) 

For the sake of simplicity we have neglected the Yukawa interactions of neutrinos. Assuming 
that the right-handed neutrinos are massive (Mjv > 10^'^ GeV) our formulas are valid below 
Mjv to a very good approximations even if such interactions are included. 
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III. RGES AND ANALYTIC FORMULAS FOR RUNNING MASSES 



Using the definition (|2.1j ) and ( p.4| )-( pTT2 ), we obtain the RGEs for the mass matrices in 



the broken phases of SM, 2HDM, or MSSM in the following form, 



167r 



dt 
^dME 



dt 

where the coefficients in the RHS are defined for the three cases. 
SM, 2HDM 



a = (?,o,8 



J2Ci'gf + cVEY^^ME (3.1) 



5' 

a= f-J>o,8 

Cr= (|,0,0) (3.2) 



MSSM 



49 9 16 

20'4'y 
19 9 16 

60'4'y 



.20 4 

=(3,1,3) (3.3) 



SM 

'3 3 3 



2HDM 



(a.M)^(|-y) (3.4) 



(«.M)^(|4,|) (3.5) 
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Defining the diagonal mass matrices Mp, the diagonal Yukawa matrices {Yp) and the CKM 
matrix {V) through biunitary transformations Lp and Rp on the left(right)-handed fermion 
Fl{Fr) with F = U,D,E 

Mp = L^pMpRp 
Yp = L^pYpRp 
Ml = L^pMpM^Lp 



^ F 



L^pYpYpLp 



V = lILd 



(3.6) 



and following the procedures outlined in 0,^] we obtain 



dt 



dM 



D 



dt 

dMl 
dt 



mIj,lIlu 



Ml^L^pLp 



+ 



+ 



1 

167r2 
1 

16^ 



-2 CigjMl + 2~aY^Ml + h [VYIV^MI + MIVYIV^ 

i 

-2 C[g]Ml + 2~aYlMl + h (v^Y^VMl + MlV^Y^V 



-2YC:'9^M'E + 2cYiMl 



(3.7) 



where Lp — . 

We point out that in the corresponding RGEs for Yukawa couplings given by eq. (2.13) in 
Ref. 0, the terms -2 C^g^Y^I (levr^), -2 CfgfY^I (IGtt^) and -2 Q^^f^l/ (levr^) 
are missing from the R.H.S. 

The diagonal elements of L^pLp{F = U,D,E) are made to vanish in the usual manner 



by diagonal phase multiplication. The nondiagonal elements of both sides of (|3.7|) give the 
same RGEs for CKM matrix elements as before which on integration yields p8| , |29[| , 



|K,/3(mt)| e-i^(-^*('^)+^''(^)), a(3 = ub,cb,tb,ts 



(3.8) 



otherwise 



IK/3 ("it) I ) 

Taking diagonal elements of both sides of ( p.7| ) and using dominance of Yukawa couplings 
of the third generation over the first two, except the charm quark, we obtain RGEs for the 
mass eigen values of quarks and leptons. 
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IGtt 



dt 

.drrir 



-Y^Qgf + byllV. 



ub 



dt 

2dmt 



-J2ag^ + ~ayl + byl\V, 



cb\ 



dt 
,dmi 



-J2ag^ + ~ay^ + by',\Vt,\ 



rrir 



rrit 



IGvr 



dt 
2dmb 



Y.C'd + hyl\V,, 



mj,j = d,s 



dt 
drrii 



rrib 



IGtt 



dt 



f^j,] = e,/i 



dt 



(3.9) 



Integrating ( p.9| ) and using the corresponding low-energy values, the new analytic formulas 
are obtained in the small mixing limit as 

m„(/i) = m„(l GeV)?7^^5^^ 
mc(/i) = mcirricj'n'^B-^e''^'' 
mt{iJ.) = mt{mt)B-'e'''^+''^ 
mi(fi) = mj(l GeV)r7,^^-B^"'^, i = d,s 
rribiiJ.) = mb{mb)%'B^'e^'^+''^ 
mi(fi) = mj(l GeV)ri^^ , i = e,/i 
nirin) = mr{m^)ri~^B~^e''^^ 



(3.10) 



where 



B. 



n 
n 
n 



ai{mt) 



2b,- 



C' 
i_ 

2b,- 



C" 

%_ 

2b, 



(3.11) 
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The ratio rif{f = u, d, c, s, b, e, /x, r) appearing in ( p.lO| ) is the QCD-QED rescahng factor for 
fermion mass ruf. 

Comparison with earher and recent derivation of analytic formulas 0,0] shows several 



differences. Whereas the top-quark Yukawa coupling integral in (|3.12| ) has been predicted 



to affect the running of and mc(yu) our formulas predict no such dependence. Simi- 

larly, whereas the 6-quark and the r-lepton Yukawa coupling integrals have been predicted 
to affect the running charged- lepton masses me(/u) and m^(yu), our formulas predict no such 
contributions. Our formulas predict that in all the three cases, SM, 2HDM, or MSSM, the 
third generation Yukawa couplings do not affect the running masses of the first two gener- 
ations in the small mixing limit. Also for the third generation running masses, the Yukawa 
coupling integrals occur in the exponential factors with different coefficients as compared 
to |]30| , |3l| . The dependence on gauge couplings are also noted to be quite different in our 
analytic formulas. Whereas earlier derivation ||30| predicted the occurrence of the exponents 



Cf/26j, Cf/2bi, Cf/26j on the RHS of our formulas predict the corresponding ex- 

ponents to be Ci/2bi, C-/2bi, Cf/26j respectively. Thus, the new analytic formulas derived 
here at one-loop level predict substantially new dependence on gauge and Yukawa couplings 
for the running masses. Our formulas for the case of MSSM are the same as those obtained 
in []13| where the SUSY scale was taken as Ms = rrit. As derived in for the MSSM, the 



formula for running tan/3(/i) has also the same form also in the 2HDM at one-loop level. 

tan/9(/i) = tan P{mt) exp[—3It{fJ.) + 3Ib{fJ.) + IrifJ')] (3.13) 

In contrast, when running of VEVs are ignored, apart from the mass predictions being 
different, tan/5 is also predicted to be the same for all values of /i > Mz in 2HDM or MSSM 
30,|3ll. 
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IV. FORMULAS IN MSSM FOR Ms > Mz 



In the MSSM the natural SUSY scale {Ms) could be very different from the weak scale 
with Ms ~ O (TeV), whereas Ms ^ 1 TeV has gauge hierarchy problem. As our new 
contribution in MSSM in this paper, compared to we present new analytic formulas for 
all charged fermion masses for any SUSY scale Ms > Mz by running them from mt — Ms 
as in SM and then from Ms — /i as in MSSM. 

m„(/i) = m„(l GeY)^]-^Gu{^^) (4.1) 

'3 



m,(/i) = m,(m,)r/-^G„(/i) exp (^-Ic{Ms) + 3/c(Ai)) (4.2) 

mt{^l) = mt{mt)GM exp {^It{Ms) - \h{Ms) + + h{^^)) (4.3) 

mi{iJi) = mi{l GeV)ril'^Gd{li),i = d,s (4.4) 

mf,(/i) = mb{mb)v;'Gdifi) exp (^hb{Ms) - ^ItiMs) + Shifi) + /t(/i)) (4.5) 

mi{jj) = mj(l GeV)r7~^G'e(/i), i = e, /i (4.6) 

rririfi) = m^{m^)7]~^Ge{^) exp {^Ir{Ms) + 3/^(;u)^ (4.7) 



where 

GM 
Gdifi) 



-2 4 -49 -9 

ai(M5)V' as{Ms)Y ai(/i) / a2(/i) \ « / a^il^) 



aiimt) ) \a3{mt) J \ai{Ms) J \a2iMs) J yasiM; 



s 



' ai{Ms) \ « f a3{Ms) \ ^ ( aM \ ( aM_\ ~' ( aM 
^aiimt) \asimt) [aiiMs) ^2(^5)/ [asiMs] 



ai 



and If {Ms) is defined through (|3.12| ) with /i = M5. Running of the elements of the CKM 
matrix in the MSSM are modified by the following formulas 

|K/3(mt)| ei(^*(*^s)+4{Afs))e-(^*(^)+^''>(A')), = ub,cb,tb,ts 
KM\ = { (4.10) 
iVapirrit)] , otherwise 
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The one-loop formula for tan/3(yu) in ( |3.13| ) is also modified 



tan/3(/i) = tan (3{Ms) exp (-3/t(/i) + 3h{fi) + /.(/i)) (4.11) 

The analytic formulas ( [4.1|) -( p?TD hold good for any value of < Ms < /i. It may be 
noted that in the limit of Ms rrit, If (Ms) — > 0, //(/u) -^/(/w) and the formulas ( [4.1|) - 
( |4.11| ) reduce to those obtained in [1^ . It is to be noted that corresponding exponent in the 
expression Bu in eq. (3.6) of ref. [13] should be corrected as ^ in place of 

V. NUMERICAL PREDICTIONS AT HIGHER SCALES 

The analytic formulas given in the previous section predict masses and the CKM matrix 
elements upto one-loop level at higher scales. We have also estimated numerically the effect 
of scale dependent VEVs on predictions of the running masses at two-loop level. We solve 
RGEs for the Yukawa matrices and VEVs including two loop contributions in SM and 
MSSM p!5|- p^ , |25| - p7|] numerically and obtain the mass matrices at higher scales from the 
corresponding products of the two. For this purpose the elements of the CKM matrix at 
higher scales have been obtained by running them through one-loop RGEs given by (|3.8| ) with 
appropriate values of the coefficient c given in (|2.10|) - (|2.12|) ||28|p9| . In 2HDM we carry out all 
numerical estimations at one-loop level. We use the following inputs for the running masses 
(mj), SM gauge couplings (01,02,03), electromagnetic finestructure constant (o), electroweak 
mixing angle and the CKM matrix {V) at n = Mz which have been obtained from the 



experimental data [32 



TTiu = 2.33toil MeV,me = 677^^? MeV, 
mt = 181 ± 13 GeV,mrf = 4.69i;^|[? MeV, 
rus = QSAtll^Q MeV, nib = 3.00 ± 0.11 GeV, 
me = 0.48684727 ± 0.00000014 MeV, 

= 102.75138 ± 0.00033 MeV, 
rrir = 1.74669l°:S GeV (5.1) 
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ai{Mz) 
a2{Mz) 



0.016829 ±0.000017 



0.033493 



-0.000042 
^-0.000038 

0.118 ±0.003 
128.896 ±0.09 



V{Mz 



sin^ Ow = 0.23165 ± 0.000024 



0.9757, 0.2205, 



-0.2203 - O.OOOle*^ 0.9747, 



0.0030e-^'5 
0.0373 



0.0082 - 0.0029e^ 



-0.0364 - 0.0007e^^ 0.9993 



(5.2) 



(5.3) 



For the sake of convince we have used 5 = 7r/2 as in The choice of same inputs enables 
us to compare our results on mass predictions with those obtained with scale-independent 
assumption on the VEVs in SM and MSSM |3^. We neglect mixings among charged-leptons 
and use the diagonal basis for up-quarks. 

The variations of VEVs as a function of fi are shown in Figs. |I|-^for the SM, 2HDM, and 
MSSM where the initial value of tan/5(M^) = 10 has been used for the latter two cases. In 
these and certain other Figs, we have used the variable t = In /i along X-axis where /i is in 
units of GeV. It is quite clear that in the SM as well as the other cases the running effects 
of the VEVs contribute to very significant departures from the assumed scale-independent 



values P8| , p0| -p^. Thus the predicted running masses are to be different in all the three 
cases. Since fu(/i) increases and fd(/i) decreases with increasing /x, the up-quark masses are 
expected to have decreasing effects whereas the down-quark and charged-lepton masses are 
expected to have increasing effects at higher scales in MSSM and 2HDM. But in the SM 
all the masses are expected to have decreasing effects due to decreasing value of f (/i). In 
fact these features are clearly exhibited in all numerical values of mass predictions carried 
out in this investigation. It is to be noted that almost all fermion masses, except the top- 
quark, 6-quark and the r-lepton near the perturbative limits, decrease at higher scales due 
to decrease in corresponding Yukawa couplings. But the effect of running VEVs contribute 
to additional decreasing or increasing factors in the respective cases. 
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The predictions of all the charged fermion masses as a function of t = In /i are shown 
in Fig. I with Ms = Mz and in Fig. | with Ms = I TeV in the case of MSSM using 
ta.n j3{Mz) = 10. The corresponding predictions in 2HDM and SM are shown in Figs. |^ 
^. In Fig. we display the comparison of mass predictions as functions of t = In/x with 
and without running VEVs in MSSM assuming Ms = Mz and tan f3 = 10. Although 
the differences in the two types of predictions are clearly distinguishable, they are quite 
prominent in the up-quark sectors. While the new contributions are seen to be significant 
for the down-quarks and charged-leptons at higher scales with > 10^ GeV, in the case 
of up-quarks the contributions are found to be important starting from fi = O (TeV). As 
compared to the scale-independent assumption |3^, our predictions are clearly smaller for 



the up-quarks and larger for the down-quarks and charged-leptons as indicated by solid-line 
curves in Fig. |^. With the input values for nit and rrif, in ( |5.1j ), the lowest allowed value 
of tan/?(Ms) is determined by observing the perturbative limit for the top-quark Yukawa 
coupling at the GUT scale, uI^Mgut) / '^t^ < 1.0 and the highest allowed value of taxi [3 {Ms) 
is determined from the corresponding limit on the 6-quark Yukawa coupling. 
MSSM 

Ms = Mz: 2.3lf)i < tan (3{Ms) < 58.7^|^, (5.4) 
Ms = 1 TeV : l.7tli < tan l3{Ms) < 64.8+|^ (5.5) 

The allowed region for tan l3{Ms) as a function of Ms in MSSM is shown in Fig. ^ where 
the solid (dashed) lines are due to the central values (uncertainties) in the inputs of rrit and 
rrih. It is clear that the allowed region for tan (3 increases, although slowly, with increasing 
Ms- In the 2IIDM the allowed region for tan/3 is found to be substantially larger. 
2HDM 

1.2+°;^ < tan/5(Mz) < 68.9 ±2.7 (5.6) 

We have noted that in all the three effective theories the difference between one- and 
two-loop estimation of running masses at the highest scale (M{/) varies between 1-5%, the 
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lowest discrepancy being for the leptons and the highest being for the top-quark. But in 
MSSM and 2HDM this discrepancy increases to 10-12% for the b- and the top-quarks as the 
respective perturbative hmits are approached. 

The running VEVs in MSSM and 2HDM lead to variation of tan/?(yu) as a function of 
over its initial value at Mz- This is shown in Fig. ^ for different input values where the 
dashed (solid) line represents the case for 2HDM (MSSM). In both the theories tan/?(/i) 
decreases (increases) from its initial value when the latter crosses a critical point. This 
critical value is tan p{Mz) ^ 56 (52) in MSSM (2HDM). In Fig. we present tan p{Mu) 
at the GUT scale as a function of tan j3{Mz) for both the theories. We observe steep rise 
in the curves as the respective perturbative limits are approached in the large tan (3{Mz) 
region. 

Using the central values of mt{Mz), mh{Mz) and itIt-^Mz) from (|5.1|) , we have studied 
variation of mt{fi), mb{fi) and mr{^) for different values of /i = 10^ GeV, 10^^ GeV and 
2 X 10^^ GeV each as a function of various low-energy input values of tan (3{Mz) in MSSM 
and 2HDM. These results are presented in Figs. |TT]-|I3| for the 2IIDM (dashed lines) and for 
the MSSM (solid lines) with Ms = Mz- It is clear that the perturbatively allowed range of 
tan (3 decreases with increasing /i both for MSSM and 2HDM. 

We have examined simultaneous variation of mt{fi) as a function of fi and tan/?(/i) which 
is displayed in the three dimensional plot of Fig. |14| for the input value of mt{Mz) = 181 
GeV and tan (3{Ms) = 2 — 58, where Ms = 1 TeV. Using the top-quark mass at /i = Mz, 
we have calculated mt{^) and tan/5(/x) at every /i between Ms — Mjj for input value of 
tan j3{Ms) = 2 — 59. The results are displayed in the three dimensional plot. The variations 
of the running mass predictions at the GUT scale {Mu = 2 x 10^^ GeV) as a function of 
the SUSY scale {Ms = Mz — 10^ GeV) are shown in Figs. p!5| - |T6| for the third generation 



fermions using various input values of tan j3{Ms). We find that the top-quark mass at the 
GUT scale at first decreases sharply in the smaller and larger tan (3 regions as Ms increases 
and then remains almost constant for Ms = 3 x 10^ — 10^ GeV. Similarly the predicted b 
and r masses decrease with increasing Ms although fall off is slower in the case of r in the 
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large tan /3 region. 

Numerical values of predictions of the running masses are presented in Table. | for the SM 
at three different scales fi = 10^ GeV, 10^"^ GeV and 2x 10^^ GeV. The two-loop contributions 
to the RG evolution of Yukawa couplings depends, although very weakly, upon the Higgs 
quartic couplings A which is related to the Higgs mass {M^) and VEV(f), A = Mfj/v"^. We 
have used the two- loop RGEs for A(/i) for the SM |jl6| and evaluated the running masses 



and VEVs of Table. | for the input value of the Higgs mass Mh = 250 GeV. Changing the 



Higgs mass in the allowed range of Mh = 220 — 260 GeV ^2| does not change the results 
of Table. | significantly. The uncertainties in the quantities are due to those in the running 
masses at /i = Mz- The mass matrices for Mfe(yu) and M„(yu) are modified by the factors 
v{^) /v{Mz) where v{Mz) ~ 174 GeV and the CKM matrices at higher scales remain the 
|3^ . The computed values of masses are found to be less when compared to those 



same as m 



obtained with scale- independent assumption ||32|. This is clearly understood as the running 
VEV in the SM decreases with increasing fi. For example in the SM at = 2 x 10^^ GeV 



our predictions are m^, m^) = (0.83 MeV, 242.6 MeV, 75.4 GeV) as compared to 



(m„,mc,mt) = (0.94 MeV, 272 MeV, 84 GeV) where the running effect on the VEV has 
been ignored. In Tables [ni-[m| numerical values of masses, VEVs and tan (3 are given at the 
same three scales for the MSSM and 2HDM with tan f3{Mz) = 10 and 55 in each of the two 
theories. As emphasized in this paper our high scale estimations predict quite significantly 
different values for the running masses especially in the up-quark sector. Although the 
CKM matrices at high scales remain the same as in scale-independent assumptions on the 
VEVs, the up-quark mass matrices are modified by the factor Vu{fJ')/vu{Ms), but the down- 
quark and charged-lepton mass matrices are modified by the factor Vd{fi)/vd{Ms). In the 
MSSM with Ms = Mz and tan (3{Mz) = 10 including running effect of the VEVs, the 
GUT scale predictions are {mujincmt) = (0.70 MeV, 200 MeV, 73.5 GeV) as compared 
to |2[ {mu,mc,mt) = (1.04 MeV, 302 MeV, 129 GeV). But by increasing the SUSY scale 
to Ms = 1 TeV and in the large tan /3-region we find substantial decrease in the predicted 
values of the top-quark mass at the GUT scale leading to {rriu^rric^mt) = (0.72 MeV, 210 

16 



MeV, 95.1 GeV). This is understood by noting that tan/3 ^ 55 is closer to the perturbative 
hmit for which the top-quark Yukawa couphng is larger. Similarly from Table. ^ we note 
nearly a 20% increase in the mb{Mu) with the increase of tan j3 from 10 — 55. Similar effects 
are also noted in 2HDM as can be seen from Table. ^ where rrit decreases by nearly 7% as 
tan (3 increases from 10 — 55. For larger effect the increase has to be larger in tan (3 since the 
perturbative limit in this case is closer to tan/? ^ 69 as compared to the MSSM case where 
the limiting value is tan/3 ^ 59. 



VI. SUMMARY AND CONCLUSION 

In the frequently exploited renormalisation scheme in gauge theories, the Yukawa cou- 
plings and VEVs in the SM, 2HDM and MSSM run separately The effect of scale 



dependence of the VEVs has been ignored while deriving analytic formulas |28,30,31| and 



obtaining numerical predictions at higher scales for the running masses of fermions [32 



but appropriately taken into account more recently [O]. In this paper we have derived 



new analytic formulas in the SM and 2HDM and generalised the formulas of [|1^] for any 
supersymmetry scale {Ms > Mz)- The new formulas exhibit substantially different func- 
tional dependence on gauge and Yukawa couplings in all the three effective theories. In 
particular, the running masses of the first two generations are found to be independent of 
Yukawa couplings of the third generations in the small mixing limit. Numerical predictions 
at two-loop level shows that all the running masses in the SM and only the up-quark masses 
in the MSSM and 2HDM decrease at high scales when compared with the predictions tak- 
ing scale-independent VEVs. But in the case of MSSM and 2IIDM the down-quark and 
the charged-lepton masses increase over the corresponding predictions obtained with scale- 
independent assumptions on the VEVs. Compared to the MSSM the perturbatively allowed 
region of tan (3 is larger in 2IIDM. In MSSM the allowed region shows a slow increase with 
the SUSY scale. We have also made predictions of the running masses at the GUT scale as 
a function of supersymmetry scale exhibiting new behaviours. We suggest that these new 
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analytic formulas and improved estimations on the running masses and tan (3 at high scales 
be used as inputs to test models proposing unified explanations of quark-lepton masses. 
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FIGURES 
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FIG. 1. Variation of running VEVs in the SM, 2HDM and MSSM as a function of fi{t = In/x) 
showing substantial deviation from the scale-independent assumption. 
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FIG. 2. Variation of running VEVs at higher scales in MSSM and 2HDM as a function of 
IJ,{t = In f^) showing substantial deviation from the scale-independent assumption. 
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FIG. 3. Predictions of running masses at higher scales as a function of (t = Irifi) in MSSM 
with SUSY scale Ms = Mz using the input parameters given in (^.1|)- (|5.3| ) and tan (3{Ms) = 10. 
The dashed lines are due to uncertainties in the input parameters. 
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FIG. 5. Predictions of running masses at higher scales in the 2HDM using the input parameters 
given in Q-Q and tan p{Mz) = 10. 
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FIG. 6. Predictions of running masses at higher scales in SM with the input parameters given 
in Q-Q and Mh = 250 GeV. 
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FIG. 7. Comparison of running mass predictions in the MSSM (solid lines) with those obtained 
from scale-independent assumptions (dashed lines) on the VEVs. The SUSY scale has been taken 
to be Mz- 
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FIG. 8. Perturbatively allowed region for tan P{Ms) as a function of SUSY scale Ms- The 

lower (upper) limits are due to top-quark (6-quark) Yukawa coupling. The dashed lines are due to 

uncertainties in the respective input masses. 
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FIG. 9. Variation of tan as a function of fi (t = In fi) for different input values of tan (3{Mz) 
in MSSM (solid lines) and 2HDM (dashed lines). The inputs for different curves starting from the 
bottom most Une is tan/3(Mz) =2, 5, 10, 20, 30, 40, 50 and 58. In the MSSM the SUSY scale has 
been taken as Mz 
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FIG. 11. Prediction of mt{p) at higher scales, n = 10^ GeV, 10^^ GeV and 2 x lO^^ GeV as a 



function of tan l3{Mz) in MSSM with Ms = Mz (sohd hnes) and 2HDM (dashed hnes). 
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FIG. 12. Predicton of mbip,) at higher scales, fi = 10^ GcV, 10^^ GeV and 2 x lO^^ GeV as a 
function of tan l3{Mz) in MSSM with Ms = Mz (sohd hnes) and 2HDM (dashed lines). 
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FIG. 13. Predicton of m^(/i) as at higher scales, fj. = 10^ GeV, 10^^ GcV and 2 x 10^^ GeV as 



a function of tsm l3{Mz) in MSSM with Ms = Mz (solid lined) and 2HDM (dashed lines). 
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FIG. 14. Prediction of top-quark mass mt{fi) at higher scales {/j, > Mz) as a function of 
{t = In/i) and tan/3(/u) in MSSM with Mg = 1 TeV. The values of tan/3(/i) at very ji has been 
obtained through solutions of the corresponding RGE using iw.f3{Mz) = 2 — 58 as inputs. 
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FIG. 15. Variation of top-quark mass prediction at the GUT scale as a function of SUSY scale 
Ms = Mz- 10^ GeV and various values of tan/3(M5)=3 (solid line), 10 (large-dashed line), 50 
(small-dashed line), 55 (dotted line). 
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FIG. 16. Same as Fig. 15 but for 6-quark and r-lepton mass predictions. The input values of 
tan P{Ms) for four different curves in each case are tan P{Ms) = 3, 10, 50, and 55 in the increasing 
order of masses. 
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TABLES 

TABLE L Running mass and VEV predictions at higher scales in the nonSUSY standard model 
for the input values of the Higgs mass Mh = 250 GeV and other parameters given in (5.1)-(5.3). 
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TABLE II. Predictions of running masses, VEVs and tan (3 at higher scales n = 10^ GeV, 10 
GeV and 2 x 10^^ GeV in MSSM with SUSY scale Ms = 1 TeV, using two-loop RG equations. 
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tan P{Ms) = 55 


fi = 10^ (GeV) 
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TABLE III. Predictions of running masses, VEVs and tan P in 2HDM at higher scales using 



one-loop RG (;quMlious. 
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tan P{Ms) = 55 


fi = 10^ (GeV) 


^ = iQi^ (GeV) 


fx = 2x 10^^ (GeV) 


rriu (MeV) 


1.2021+°-2309 

" — 0.2417 


0.9908+^-i^?,^ 

w.^iywu_o 2002 


0.8749+^-}^Si 

vj.u 1 ^»'_o.l772 


rric (MeV) 


349.2889+'H-?Eg^ 

34.5782 


287.9066+ol-oo^^ 

28.8338 


254.2223+^tmtl 

" 25.6031 


rrif (GeV) 


106 6700+Jn nlQ? 


92 1000+1^1?^^ 


oq qq-i y+13.2279 
oo.vo± 1 _io.3226 


rriri fMeV) 


■^•^"^ 1 —0.3660 


2 0430+n^nRS 
■^■"^"-"-'—0.3069 


1 8204+°-^°^° 

j..u.i.vjT:_o 2743 


rus (MeV) 


48.8888+E-o?^Q 


40.6898+«n5Qfi 

6.0498 


36.2584+!-5noo 

5.4099 


till, fCioVl 


1 1 -'-■'-(). l()<)2 


1 5392+"- j272 

--■ — ().1()<)2 


-■ -'- —(). li()2 


me (MeV) 


0.4662+0-OOOi 

" 0.0001 


0.4529+°-ZI 

0.0001 


0.4407+°nm°^ 
" " 1 0.0001 


m„ fMeV) 


98.4302+°-°°^^ 


'^"•'-'•^""—0.0094 


93.0536+^-^jf^ 

"^"■"""'-'—0.0136 


rur (GeV) 


l-7659ingi 


1 777t;+0.0073 
1. 1 ' 1 <J— 0.0060 


1 ysc^l +0-0136 
1. ( ool_o 0107 


tan P 


54.9963;ii|^ 


55.7094;2;4787 


^8'^1 -3.2730 

00.0001^;^ gggg 


vu (GeV) 


1 71 78-7-8644 
100. <±(0^g 5265 


152.0846;^2gl^^i 


150.4478;!!^^^^ 


Vd (GeV) 


QQ1/I-0 0649 
Z. 5014^0.0578 


2 72qq-°i°^2 

z. /zyy^O.0892 


2.6588;°:}^°^ 
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